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^ ■ Abstract 
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CN I General arguments in favor of the necessity of a wave packet description 

^ • of neutrino oscillations are presented, drawing from analogies with other 

Q ■ wave phenomena. We present a wave packet description of neutrino 
oscillations in stationary beams using the density matrix formalism. 

Cf) ', Recent claims of the necessity of an equal energy of different massive 

_ neutrinos are refuted. 
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1 Introduction 



The physics of massive and mixed neutrinos is one of the hot topics in today's research 
in high energy physics, following the discovery of oscillations of atmospheric and solar 
neutrinos (see the reviews in Refs. Ii3,|4il5lli3 and references therein, as well as the 
references in |4fll|). 

Neutrino oscillations have been proposed in the late 50 's by Pontecorvo 

@,E3- The 

oscillations are generated by the interference of different massive neutrinos, which are 
produced and detected coherently because of their very small mass difference. 

In 1962 Maki, Nakagawa and Sakata js^ considered for the first time a model with 
mixing of different neutrino flavors. In 1967 Pontecorvo proposed the possibility of solar 
neutrino oscillations j^, before the discovery in 1968 of the solar neutrino problem in 
the Homestake experiment ji^. In 1969 Pontecorvo and Gribov considered 
oscillations as a possible explanation of the solar neutrino problem 4^ 1 . 

The theory of neutrino oscillations in the plane-wave approximation was developed 
in the middle 70's by Eliezer and Swift jl^, Fritzsch and Minkowski [2^, Bilenky and 
Pontecorvo ^3], and beautifully reviewed by Bilenky and Pontecorvo in Ref. [l5l |. 

In 1976 Nussinov for the first time considered the wave packet nature of prop- 
agating neutrinos and inferred the existence of a coherence length, beyond which the 
interference of different massive neutrinos is not observable. This is due to the different 
group velocities of different massive neutrinos, that causes a separation of their wave 
packets. In 1996 Kiers, Nussinov and Weiss first pointed out the importance of the 
detection process for the coherence of neutrino oscillations and discussed some implica- 
tions for the wave packet approach. 

In 1981 Kayser [s^l presented the first detailed discussion of the quantum mechanical 
problems of neutrino oscillations, pointing out the necessity of a wave packet treatment. 

Wave packet models of neutrino oscillations have been later developed in the frame- 
work of quantum mechanic s IsTilsii E3 . Isi , 2^, 23l.l83l| and in the framework of quantum 
field theory 3^ 3^ 54, 3, id l32l| fsee also Ref. j55l] and the reviews in Refs. 74, ll|). 

In spite of the well-known fact that in quantum theory localized particles are described 
by wave packets and in spite of the success of the wave packet treatment of neutrino oscil- 
lations, some authors have presented arguments against such approach. In Section El we 
briefly review well-known evidences of the wave packet nature of light, which, by analogy, 
imply a wave packet nature of neutrinos. In Sections El and El we present objections to 
arguments against a wave packet description of neutrinos. In Section HI we present also 
a wave packet derivation of neutrino oscillations in the density matrix formalism and we 
discuss some implications for the case of stationary beams. 

Another debated problem in the theory of neutrino oscillations is the determination of 
the energies of massive neutrinos. In SectionlHlwe present objections to recently proposed 
arguments in favor of the equality of the energies of different massive neutrinos. 



2 No source of waves vibrates indefinitely 

The title of this section is taken from the beginning sentence of Section 11.11 of Ref. 
which concludes with the following paragraph: 
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In light sources, the radiating atoms emit wave trains of finite length. Usually, 
because of collisions or damping arising from other causes, these packets are 
very short. According to the theorem mentioned above^, the consequence is 
that the spectrum lines will not be very narrow but will have an appreciable 
width AA. A measurement of this width will yield the effective "lifetime" 
of the electromagnetic oscillators in the atoms and the average length of the 
wave packets. A low-pressure discharge through the vapor of mercury contain- 
ing the single isotope ^^^Hg yields very sharp spectral lines, of width about 
0.005 A. Taking the wavelength of one of the brightest lines, 5461 A, we may 
estimate that there are roughly 10^ waves in a packet and that the packets 
themselves are some 50 cm long. 

The broadening of optical lines due to the finite lifetime of atomic transitions is known 
as natural linewidth. 

It is interesting to note that the broadening of optical lines was well known to exper- 
imental physicists in the nineteen century and explained by classical models before the 
advent of quantum theory (see Ref. references therein). In 1895 Michelson 6^ 



listed among the hypotheses formulated before that time to account for line broadening 

3. The exponential diminution in amplitude of the vibrations due to commu- 
nication of energy to the surrounding medium or to other causes. 

As explained in [l7| . 

We consider an emitting atom which we shall proceed to remove to infinity 
and reduce the "temperature" to the point where, classically at least, no 
translational motion exists. Now from the classical picture of a vibrating 
electron or the simple picture of a pair of energy levels between which our 
radiation transition takes place, we should expect these conditions to yield a 
spectral line of a single frequency. We, of course, do not obtain this result, 
but, rather, we obtain the familiar natural line shape which is attributable to 
Michelson's Cause 3. 

A classical derivation of the natural linewidth can be found in Section 17.7 of Ref. 0|. Let 
us emphasize that the natural linewidth of atomic lines has been observed experimentally 
(see Section 21.4 of Ref. 0)! 

Another important cause of line broadening known in 1895 was 

4. The change in wavelength due to the Doppler effect of the component of 
the velocity of the vibrating atom in the line of sight. 

The Doppler line broadening due to the thermal motion of atoms in a medium was 
calculated by Rayleigh in 1889 ji^. This important effect, that must be aways taken into 
account in calculating the spectral shape of monochromatic beams, does not concern us 
here, because it does not generate a coherent broadening. It is simply due to the different 
motion of different atoms, whose radiation is incoherent. 

To these (and other) causes of line broadening Michelson added in 1895 j6o| 

^ "The largest the number N of waves in the group, the smaUer the spread AA, and in fact theory 
shows that AA/Ao is approximately equal to 1/iV." (Section 11.11 of Ref. [i^.') 
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5. The limitation of the number of regular vibrations by more or less abrupt 
changes of phase amplitude or plane of vibration caused by collisions. 



This important coherent effect has been called with several names, among which collision 
broadening, pressure broadening and interruption broadening, it has been studied in depth 
by many authors (see, for example, Refs. jlTllHill) and it has been observed experimentally 
(see Section 21.5 of Ref. 0)! 

In quantum theory, the fact that "no source of waves vibrates indefinitely" implies 
that all particles are produced as wave packets, whose size is determined by the finite 
lifetime of the parent particle, or by its finite mean free path if the production process 
occurs in a medium. Since also no wave detector vibrates indefinitely, it is clear that all 
particles are also detected as wave packets. 

Usually, at least in high energy physics, the wave packet character of particles is 
not important and they can be well approximated by plane waves. The phenomenon 
of neutrino oscillations is an exception, as shown in 1981 by Kayser [HH, because the 
localization of the source and detector requires a wave packet treatment. 



The wave packet approach to neutrino oscillations implemented in Refs. 87L 38L l24l . 
in the framework of quantum mechanics and in Refs. [36|, |39|, |5J, |18|, |10|, |3^ in 



the framework of quantum field theory allows to derive the neutrino oscillation probability 
from first principles in a consistent framework. 

In the light of these considerations, it seems rather surprising that some physicists 
do not agree with a wave packet description of neutrinos and put forward arguments in 
favor of a unique value of neutrino energy. A tentative explanation of this approach stems 
from the common education and practice in physics of working with energy eigenstates, 
forgetting that they are only approximations of the states of physical systems in the real 
world^. In the following Sections we present a critical discussion of recently advocated 
approaches to neutrino oscillations in which neutrinos are not described by wave packets 
with some spread in momentum and energy. 



3 Bound states 

Bound states of atomic electrons are usually calculated with different degrees of ap- 
proximation taking into account the static potential generated by the nucleus and other 
electrons and neglecting the coupling of the electron with the electromagnetic radiation 
field. In this approximation, bound states are stationary and have definite energy. In 
spite of the obvious fact that in such approximation electrons cannot jump from one 
bound state to another, the energy difference between two stationary bound states is 
sometimes associated without any further explanation with the energy of the photon 
emitted or absorbed in the transition of an electron from one bound state to the other. 

However, in order to describe the transitions between bound states it is necessary to 
take into account the coupling of the electron with the time-dependent electromagnetic 
radiation field. Since the bound states are not eigenstates of the full Hamiltonian, in 
reality they are not stationary and they do not have definite energy. The lack of a 
definite energy is a necessary requirement for the existence of transitions: since the 

would like to thank I. Pesando for an illuminating discussion about this point. 
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Schrodinger equation implies that the time evolution of energy eigenstates evolve by a 
phase, orthogonality between different energy eigenstates is conserved in time, excluding 
transitions. 

The fact that unstable systems do not have a definite energy is well-known to high 
energy physicists from their experience with resonances (hadronic excited states) and 
highly unstable particles {Z and W bosons), which do not have a definite mass (energy 
in the rest frame). Their mass width is inversely proportional to the decay rate. 

From these considerations it follows that the quantum field theoretical model of neu- 
trino oscillations in the process 

n p + e~ + I'e ""^^^^ i^e + e~ — > z/g + e~ (3.1) 



presented in Ref. j4^, in which a neutrino is assumed to be produced by a nuclear 
/3-decay between two stationary bound states with definite energy and detected through 
scattering with an electron in a stationary atomic state with definite energy, is unrealistic, 
not to mention the description of the final-state electrons by plane waves, which is in 
contradiction with the fact that the electrons interact with the surrounding medium and 
with the necessity to observe the final lepton in the detection process, as already noted in 
[3^ . Of course, once the matters of principle are clear, one can consider the assumptions 



in Ref. as approximations acceptable in some cases. One must also notice that 



Ref. 4^ is very interesting for the useful theorem proved in the appendix, which has 
been used in the quantum field theoretical calculation of neutrino oscillations in the wave 
packet approach in Refs. 3^, 10^ . 



Developing the technique proposed in Ref. j4^, the authors of Refs. j45|,|43[ considered 
neutrino oscillations in the process 



^ e + Ue + > + P ^ n + e (3.2) 

with the muon correctly described by a wave packet, but the proton was still described 
by an unrealistic stationary bound state with definite energy. Furthermore, all final 
particles were described by plane waves, in contradiction with the fact that the positrons 
and neutron interact with the surrounding medium and the necessity to observe the final 
positron in the detection process in order to measure oscillations. Localized interactions 
of a particle with a detector obviously imply that the particle cannot be described by an 
unlocalized plane wave and a wave packet description is necessary. 



4 Stationary beams 

The author of the very interesting Ref. [7^ wrote: 

. . . many wavepacket discussions for the coherence properties of particle beams 
are unnecessary since they deal with stationary sources; and when the problem 
is stationary, essentially all information is in the energy spectrum. 

Indeed, as proved in Ref. jH^, 

. . . under very general assumptions it is not possible to distinguish experimen- 
tally neutrinos produced in some region of space as wave packets from those 
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produced in the same region of space as plane waves with the same energy 
distribution. 

It is not clear if the purpose of Ref. j^] was to show that neutrinos are not described 
by wave packets, but it is certain that it has been interpreted in this way by some 
physicists. However, Ref. [7^ certainly does not provide any help for the derivation of 
neutrino oscillations in the framework of quantum field theory. Even in the framework of 
quantum mechanics Ref. jlO] does not go beyond the well-known discussion of neutrino 
oscillations assuming equal energies for the different massive neutrino components |57i.i4d|. 
which has been shown to be unrealistic in Refs. js^ |^ and is further criticized in the 
Section 

An interpretation of Ref. [l^] as a proof that neutrinos are not described by wave 
packets seems to stem from a confusion between microscopic and macroscopic stationarity, 
as already remarked in 3, 13 • Since most neutrino sources emit neutrinos at a constant 



rate over macroscopic intervals of time, they can be considered macroscopically stationary. 
However, the microscopic processes of neutrino production (nuclear or particle decay, etc.) 
are certainly not stationary and the single neutrino does not know if it will be part of a 
stationary beam or not. Because of the localization in space and time of the production 
and detection processes, a description of neutrino oscillations in terms of wave packets is 
inescapable. 

Of course, after the matters of principle are settled, neutrino oscillations are derived 
from first principles in a consistent theoretical framework, and an estimation of the co- 
herence length has shown that it is much longer than the source-detector distance in a 
neutrino oscillation experiment, the analysis of the experimental data can be performed 
using the standard oscillation probability obtained in the plane wave approximation. This 
turns out to be the case in all present-day neutrino oscillation experiments. 

In the following part of this section we present a wave packet derivation of neutrino 
oscillations in the density matrix formalism, which is suitable for the description of a 
stationary beam, as emphasized in Ref. [l^l- Our calculation follows a method similar to 
the one presented in Ref. [s^. 

In a quantum-mechanical wave packet treatment of neutrino oscillations a neutrino 
produced at the origin of the space-time coordinates by a weak interaction process with 
definite fiavor a (a = e, /i, r) and propagating along the x axis is described by the state 

Kix,t)) = j2u:kMx.t)\u,), (4.1) 

k 

where U is the mixing matrix, |z/fc) is the state of a neutrino with mass rrik and ilJk{x,t) 
is its wave function. Approximating the momentum distribution of the massive neutrino 
z/fc with a gaussian. 



ipk{p) = [27rap j exp 



{P - Pk] 
4aP' 



(4.2) 



where pk is the average momentum and is the momentum uncertainty determined by 
the production process, the wave function is 

Mx,t) = I -^^.(p)e^^--^^'=(^)S (4.3) 
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with the energy Ek{p) given by^ 



Ekip) 



+ m 



k ■ 



(4.4) 



If the gaussian momentum distribution ()4.2|1 is sharply peaked around the mean momen- 
tum pk (^-e. if CFp <C El{pf;)/mk), the energy -Efc(p) can be approximated by 



Ek{p) Ek + Vk {p-Pk) 



where 



is the average energy and 



Eh 



Vk = 



Ek{pk) 
dEkip) 



pl + ml 



dp 



p=Pk 



Pk_ 

Ek 



(4.5) 
(4.6) 
(4.7) 



Under this approximation, the integration over p in Eq. ()4.3p is gaussian and leads to the 
wave packet in coordinate space 



'ipk{x,t) 



p2 



-1/4 



where 



exp 



-iEkt + ipkX - 



{x - Vkty 



2< 



(4.^ 



(4.9) 



is the spatial width of the wave packet. From Eq. ()4.8|1 one can see that Vk is the group 
velocity of the wave packet of the massive neutrino Vk- 

The pure state in Eq. ()4.1|) can be written in the form of a density matrix operator as 



Pa{x,t) = \Ua{x,t)){Ua{x,t) \ . 



(4.10) 



Using Eq. ()4.8j) . we have 



paix,t) 



k,j 



-i (Ek - Ej) t + i{pk- Pj) X 



{x — Vktf' {x — Vjty 



4ar 



4ar 



\uk){u,\.{A.ll) 



This space and time dependent density matrix operator describes neutrino oscillations 
in space and time. Although in laboratory experiments it is possible to measure neu- 
trino oscillations in time through the measurement of both the production and detection 
processes (see Ref. ^4])"^, in all existing neutrino oscillation experiments only the source- 
detector distance is known. In this case the relevant density matrix operator Pa{x) is 

■^We do not make unnecessary assumptions about the values of the energies and momenta of different 
massive neutrinos. The claim that different massive neutrinos must have the same energy has been 



confuted in Refs. |35|, |3l| and in further discussed in Section |31 

''^Let us also mention the possibility of time-dependent flavor oscillations of neutrinos in the primordial 
cosmological plasma (see Ref. and references therein). 
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given by the time average of pa{x,t). Since the integral over time is gaussian, one easily 
obtains^ 



Pa[X) 



Yl ^^kUaj exp 



—I 



Vk + Vj 



{Ek -Ej)- {pk - Pj] 



X 



2 ^2 



{Vk -Vj) X 

rP2 



{Ej, — Ej] 



p2 



(4.12) 



Since this density matrix operator is time independent, it is suitable for the description 
of a stationary beam in neutrino oscillation experiments, as noticed in Ref. [7fll |. 

In order to obtain the oscillation probability, it is convenient to simplify Eq. ()4.12p 
for the realistic case of extremely relativistic neutrinos (see the discussion in Section 2 of 
Ref. In general, the average massive neutrino energies Ek can be written as 



Eu 



(4.13) 



where E is the neutrino energy in the limit of zero mass and is a dimensionless quantity 
that depends from the characteristics of the production process^. From Eq. ()4.4j) . the 
corresponding momentum in the relativistic approximation is 



Pk 



E 



(4.14) 



The generality of these relations can be understood by noting that, since the neutrino 
mass rrik enters quadratically in the energy-momentum dispersion relation ()4.4|1 . its the 
first order contribution to the energy and momentum must be proportional to ml and 
inversely proportional to E, which is the only available quantity with dimension of energy. 
Using Eqs. ()4.13|1 . ()4.14|1 and the relativistic approximation 



Vk 



1 



we obtain 



(4.15) 



AmljX f AmljX 



2E 



\4V2E^al 



with Amlj = ml 



m^ 



\^f2Eal 



(4.16) 



'''The density matrix operator in Eq. 14.12|l is normalized by Pa{x) = I, where / is the identity 
operator. 

^ The quantity cannot be calculated in a quantum mechanical framework, but its value can be 
estimated from energy-momentum conservation in the production process. The calculation of the value 
of requires a quantum field theoretical treatment (see Ref. [s^)- 
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In analogy with the production process, we describe the process of detection of a 
neutrino with flavor (3 at the coordinate x = L with the operator 



X exp 



2E V 4v^E2(tD 



(4.17) 



where is a dimensionless quantity that depends from the characteristics of the detection 
process (see footnote E)) , <J^ is the momentum uncertainty of the detection process and 



The probability of z/^ — >■ transitions is given by 



P,^^,,{L) = Tr ^p^{x)Op{x - L) 
which leads to 

Puc.^up{L) = U^kU^jU/3kU*^j exp 



I AxY{vk\Pc.{x)dp{x-L)\vk) , (4.18) 



-27Ti- 



T osc 



L 

T coh 



2^2 



T OSC 

^kj 



, (4.19) 



with the oscillation and coherence lengths 



T OSC 
^kj 



AttE 
Ami- ' 



and 



t _ iV2E\ 



^2 p2 I n2 

(^x = + ^x > 



f2^p2 , f2 d2 



(4.20) 
(4.21) 



(4.22) 
(4.23) 



This result was already obtained'' in a quantum mechanical framework in Ref. |34l. w here 
the wave packet treatment with pure states presented several years before in Ref. j37| was 
extended in order to take into account the coherence properties of the detection process, 
whose importance was first recognized in Ref. js^ . Expressions for the z/q — > z/^ transition 
probability similar to Eq. ()4.19|1 have been also obtained with wave packet treatments of 
neutrino oscillations in a quantum field theoretical framework in Refs. j3fil. l39L .18^ .lOi. l33| . 

It is important to remark that the wave packet treatment of neutrino oscillations 
confirms the standard value in Eq. ()4.20|) for the oscillation length L^^f. The coherence 
length L^°^ in Eq. ()4.21|1 is the distance beyond which the interference of the massive 
neutrinos i^k and i>j is suppressed because the separation of their wave packets when 

The only improvement in the expressions above with respect to the analogous ones in Ref. Sl IS 
the introduction of ^d, that takes into account the properties of the detection process. 
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they arrive at the detector is so large that they cannot be absorbed coherently. As 
shown by Eqs. ()4.2H1 and ()4.22|1 the coherence lengths L'j?^ are proportional to the total 
coherence size, that is dominated by the largest between the coherence sizes and 
of the production and detection process. The last term in the exponential of Eq. ()4.19|) 
implies that the interference of the massive neutrinos z/^ and uj is observable only if 
the localization of the production and detection processes is smaller than the oscillation 
length. 

The localization term is important for the distinction of neutrino oscillation exper- 
iments from experiments on the measurement of neutrino masses. As first shown by 
Kayser in Ref. [5l|, neutrino oscillations are suppressed in experiments able to measure 
the value of a neutrino mass, because the measurement of a neutrino mass implies that 
only the corresponding massive neutrino is produced or detected. 

Kayser's pl| argument goes as follows. Since a neutrino mass is measured from energy- 
momentum conservation in a process in which a neutrino is produced or detected, from the 
energy-momentum dispersion relation ()4.4|) . the uncertainty of the mass determination is 



(5m,2 = ^{2EkSEkf + {2pk5pkf - 2V2Eap , (4.24) 

where the approximation holds for realistic extremely relativistic neutrinos and cxp is the 
momentum uncertainty. If the mass of Uk is measured with an accuracy better than 
Amlj, i.e. 

I /\t77,^ I 

the neutrino Uj is not produced or detected*^ and the interference of z/^ and uj is not 
observed. 

The localization term in the oscillation probability ()4.19|1 automatically implements 
Kayser's mechanism. Indeed, it can be written as 

2nr{P^Y = e(^^]\ (4.26) 



with the momentum uncertainty 



= 4^' = 4^ + 4^ • (4.27) 

p p p 

If the condition ()4.25|1 for neutrino mass measurement is satisfied, the localization term 
()4.26|) suppresses^ the interference of z/^ and Uj. 

^The smallness of dp required to satisfy the condition H4.25I) implies that and cannot be produced 
or detected in the same process, because their momentum difference is larger than the momentum 
uncertainty, or their energy difference is larger than the energy uncertainty (for extremely relativistic 
neutrinos the energy uncertainty is practically equal to the momentum uncertainty ap). In this case only 
one of the two massive neutrinos is produced or detected. 

^This argument assumes that ^ is not too small. As discussed in Ref. jS^, in realistic cases ^ is a 
number of order one, but, at least in principle, there could be some cases in which ^ is very small, or even 
zero. In these cases Kayser's mechanism cannot be implemented in the simplified quantum mechanical 
framework presented here. However, it has been shown in Ref. that a wave packet quantum field 
theoretical derivation of the neutrino oscillation probability produces additional terms, independent from 
^ or similar quantities, which suppress oscillations when the condition H4.25|l is satisfied. I would like to 
thank M. Beuthe for a illuminating discussion about this problem. 
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In the above derivation of neutrino oscillations the wave-packet description of massive 
neutrinos is crucial in order to allow the integration over time of the time-dependent 
density matrix operator t) in Eq. ()4.11|1 . which leads to the time- independent density 
matrix operator Pa{x) in Eq. ()4.12p . Of course, as emphasized in Ref. j70|, 

A single, given density matrix can arise in different ways, especially when 
incoherence is involved. 



As already noticed in Ref. j53| the density matrix operator Pq(x) can also be gener- 
ated through an appropriate incoherent average over the energy spectrum. The same is 
obviously true for the oscillation probability ()4.19p . 

In order to illustrate this point, let us consider the simplest case of two-neutrino 
mixing. The wave packet oscillation probability ()4.19|1 becomes 



Pvc^v^iL) = - sin' 



2^9 < 1 - cos 



2E 



exp 



L 

icoh 



(4.28) 

with the oscillation and coherence lengths L"*"^ and L™^ given, respectively, by Eqs. ()4.20|1 
and 1)4.211) with Am^^ = Am^ (Arn^ is the squared-mass difference and is the mixing 
angle; see the reviews in Refs. EHElEi). An incoherent average of the probability 
in the plane- wave approximation over a gaussian energy spectrum with width (Te reads 



^incoh 



1 



■ sin^ 2^ < 1 



dE' fAm^L' 

cos ( ] exp 



2T[aE 



{ E - E') 
2al 



(4.29) 

as functions of the source- 



Figure Q shows the values of the probabilities ()4.29j) and 
detector distance L for Aw? = 2.5 x 10"^ eV^, sin^ 2't9 = 1, E = 10 GeV, and aE = (Jp = 
1 GeV. One can see that it is hard to distinguish the dashed line obtained with Eq. ()4.29|) 
from the solid line obtained with Eq. ()4.28)1 . 

This result does not mean that wave packet effects are irrelevant^°. It means that in 
practice one can calculate with reasonable approximation the decoherence of oscillations 
due to wave packet effects either taking into account the momentum spread in the calcula- 
tion of the amplitude, as done in the derivation of Eq. ()4.28p . or averaging the probability 
over the same energy spread, as done in the derivation of Eq. ()4.29|) ^^. As emphasized 
in Ref. jHl, it also means that if one does not have control on both the production and 
detection processes, one cannot know if oscillations are suppressed because of incoherent 
averaging over energy of different microscopic processes or because of decoherence due 
to the separation of wave packets. However, having good control of both the production 
and detection processes, it may be possible in the far future to reduce the causes of inco- 
herent broadening of the energy spectrum and prove experimentally that oscillations can 



^°Let us clarify here that with "wave packet effects" we mean effects due to the momentum and energy 
uncertainty of a single process associated with the interruption of the emitted wave train caused by decay 
(natural linewidth) or collisions (collision broadening) discussed in Sectional These effects are usually 
not taken into account in the calculation of neutrino energy spectra (see, for example, Ref. Q)- 

^"'^One must only be careful to notice that the energy distribution in the incoherent average of the 
probability must be normalized to one, whereas the squared modulus of the momentum distribution of 
the wave packets must be normalized to one. 
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be suppressed because of wave packet effects (in practice using an approacli similar to 
the one adopted for the measurement of the natural linewidth and collision broadening 
of atomic lines discussed in Section 12)). 

Let us finally mention that the fact that coherent and incoherent stationary beams 
are indistinguishable without a theoretical analysis is well known in optics (see Section 
7.5.8 of Ref. ^3), in neutron interferometry [io, E3| and in general stationary particle 
beams M. 



5 Energy of massive neutrinos 



In the recent Ref. |56[ the claim that different massive neutrinos must have the same 
energy li^, which was refuted in Refs. jH, EH, has been renewed. If such claim 
were correct, it would mean that, at least in a quantum mechanical framework, a wave 
packet description of neutrinos is not necessary. Indeed, in this case one can construct a 
time-independent density matrix operator from the plane wave state 

\Ux,t)) = J2u:ke''''-'''>k), (5.1) 



with pk = E"^ — m 



Paix) = \i>a{x,t)){uaix,t)\ = U*,,Uaj exp [i {pk - Pj) x] . (5.2) 

k,j 

Using this density matrix operator one can easily derive the standard oscillation proba- 
bility without need of neutrino wave packets. Therefore, in this section we consider recent 
claims in favor of an equal energy of different massive neutrinos and show that they are 
faulty. 

The author of Ref. wrote: 

. . . states with different ENERGIES ARE NEVER COHERENT in any real- 
istic experiment. States of the same energy and different momenta can be 
coherent, but may not be. 



The reason offered by the author of Ref. j56| is 

The usual detector is a nucleon, which changes its state after absorbing a 
neutrino and emitting a charged lepton, and is initially either in an energy 
eigenstate or in a statistical mixture in thermal equilibrium with its surround- 
ings. No neutrino detector has ever been prepared in a coherent mixture of 
energy eigenstates and no such detector has been proposed for future experi- 
ments. 



First, one must note that in terminology of Ref. j56|] "detector" is not a macroscopic 
device, but a nucleon hit by a neutrino. Then, the claim in Ref. jH^ is easily proved to be 
wrong using the arguments presented in Section El where we have shown that all kinds 
of waves, classical or quantistic, can be emitted or absorbed only as wave packets, i.e. 
as coherent superpositions of plane waves with different frequency (energy). If the claim 



13 



in Ref. were correct, photons, as well as other particles, could never be emitted or 
absorbed. The claim that a nucleon should be initially in an energy eigenstate has been 
shown to be incorrect in Section El Furthermore, it is a mystery why all the claims on 
energy in Ref. [s^, if they were true, do not apply also to momentum. 

Also the author of Ref. j70|] claimed that different massive neutrinos have the same 
energy: 

Energy or momentum? In mixing problems, where we have to deal with lin- 
ear combinations of particles of different mass, the question comes up as to 
whether one should deal with states of the same energy or the same momen- 
tum. Since as stated above, for stationary conditions we are to perform the 
calculation as an incoherent sum over energies, we have given the answer "en- 
ergy" . Evidently, for stationary problems it is most natural to use stationary 
wave functions ~ e~^^^. 

The meaning of this paragraph seems to be that different massive neutrinos are sup- 
posed to manage to have equal energy in order to satisfy the calculational needs of some 
theoretical physicists. Therefore, we cannot consider it a proof of anything. 

A simple way to show the absurdity of the claims that massive neutrinos must have 
equal energies is to consider the transformation properties of energy and momentum from 
one inertial frame to another, following the discussion presented in Ref. jsij ]. 

If the arguments presented in Ref. (56| were correct, in order to be produced and de- 
tected coherently massive neutrinos should necessarily have exactly equal energies in the 
inertial reference frames of both the production and detection processes (which presum- 
ably coincide with the rest frames of the initial nucleons in the production and detection 
processes, in the approach of Ref. jH^; let us skip the useless discussion of the deter- 
mination of the relevant inertial reference frame in neutrino-electron scattering). Such 
a requirement would imply that the production and detection processes have to be at 
rest in exactly the same inertial system, because, as shown in Ref. the equal-energy 
requirement is not Lorentz invariant. 

Therefore if the arguments presented in Ref. were correct it would mean that 
in practice neutrino oscillations are not observable, since it is practically impossible to 
suppress random thermal motion in the source and detector. Certainly, oscillations of 
solar and atmospheric neutrinos would be impossible. Let me remind that in this dis- 
cussion we are concerned with the oscillatory terms in the flavor transition probability 
(I4.19|l . due to the interference of different massive neutrinos, which require coherence. If 
these terms are suppressed, it is still possible to measure the distance-independent and 
energy-independent flavor- changing probability 



Y,\Ue,k?\U^kW (5.3) 



or the distance-independent flavor-changing probability calculated in Ref. j65[ if matter 
effects are im por tant j73L I6ll . Ig^ . According to the results of solar neutrino experiments 
[iqI 42, 0, El EE 1^, the solar neutrino problem is very likely due to neutrino oscilla- 
tions in matter, in which the distance-independent flavor-changing probability calculated 
in Ref. is relevant. However, recently the KamLAND experiment jl^ observed in- 
dications in favor of distance and energy dependent oscillations in vacuum of reactor 
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antineutrinos due to the same mass-squared difference responsible of solar neutrino os- 
cillations. Distance and energy dependent fiavor-changing transitions have been also ob- 
served in the Super-Kamiokande [i^, Soudan 2 and MACRO 0] atmospheric neutrino 
experiments and in the long-baseline experiment K2K j^. These experimental evidences 
show without any doubt that the arguments presented in Ref. js^ in favor of an equal 
energy for different massive neutrinos are wrong. 

Indeed, if we consider for example atmospheric neutrinos produced by decays in the 
atmosphere of highly energetic pions and muons, the rest frames of the source and de- 
tector are very different and different massive neutrinos cannot have the same energy 
in both frames. If the arguments presented in Ref. j56| were correct different massive 
neutrinos could not be either produced or detected coherently, in contradiction with the 
observed oscillations. The same reasoning applies to reactor and accelerator neutrinos, 
remembering that according to the arguments presented in Ref. jH^ coherence is possible 
only if the source and detector rest frames coincide exactly (which, by the way, is an 
absurd concept, if velocity is a real- valued quantity). 

The only rigorous and correct calculation which could be erroneously interpreted^^ as 
a proof in favor of the equal energy assumption has been presented in Ref. (see also 
Section 2.5 of Ref. 0]). The argument can be illustrated with the density matrix model 
of neutrino oscillations presented in Section |3J Instead of writing the density matrix 
operator ()4.10j) in the form 1)4.111) using the wave packets ()4.8p in which the integration 
over the momentum p in Eq. ()4.3|) has been already performed, one could write the density 
matrix operator ()4.10)1 using the wave functions ()4.3|1 as 

p^{x,t) = Y.U:,U^, I ^^,(p)^^.(p')e^(^-^')--^(^^(^)-^^(^'))V,)(z/,| . (5.4) 

As in Section I3J the time-independent density matrix operator that describes a stationary 
beam is given by the average over time of Pq,(x, t). The integral over t of exp[—i{Ek{p) — 
Ej{p'))t] yields a 5{Ek{p) — Ej{p')), which means that "interference occurs only between 
wave packet components with the same energy" (see Refs. j3fl[ FH for similar results 
in electron interferometry and kaon oscillations). This means that the time average of 
Pa{x,t) is equivalent to an appropriate incoherent average over the contributions of the 
single-energy components of the wave packets. This is the reason why coherent and 
incoherent stationary beams cannot be distinguished without a theoretical analysis, as 
discussed in Section El 

Let us emphasize, however, that this argument does not mean that different massive 
neutrinos have the same energy, or that wave packets are unnecessary. The wave packet 
nature of massive neutrinos is necessary for the existence of components with the same 
energy that produce the observable interference, and the average energies of the wave 
packets of different massive neutrino are in general different. 

^^Indeed, I made such mistake in the first version of this paper appeared in the hep-pli electronic 
archive. I am deeply indebted to M. Beuthe for his comments that helped to correct this mistake and 
my wrong criticisms of Refs. 
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6 Conclusions 



Starting from the well-known fact that "no source of waves vibrates indefinitely" , we have 
argued that neutrinos, as all other particles, are naturally described by wave packets. 
This is in agreement with the well known fact that in quantum theory localized particles 
are described by wave packets. Even when neutrinos are produced or detected through 
interactions with particles in bound states, they are described by wave packets, because 
bound states in interaction with an external field do not have a definite energy. 

Since the production and detection processes in neutrino oscillation experiments are 
localized in space-time, a consistent description of neutrino oscillations requires a wave 
packet treatment. In particular, the wave packet character of massive neutrinos is crucial 
for the derivation of neutrino oscillations in space, because the group velocity establishes 
a connection between space and time and allows the time average of the space and time 
dependent probability (or density matrix). 

Finally, we have shown that the claimed arguments in favor of an equal energy of 
massive neutrinos (which could allow a quantum mechanical derivation of neutrino os- 
cillations without wave packets) are in contradiction with well known physical laws and 
phenomena. 
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Figure 1: Oscillation probability as a function of distance L for Am^ = 2.5 x 10~^ eV^, 
sin^ 2'd = 1, and = 10 GeV. Dotted line: Unsuppressed and unaveraged oscillation 
probability. Dashed line: Oscillation probability ()4.29j) averaged over a gaussian energy 
spectrum with width cte = I GeV. Sohd line: Oscillation probability ()4.28j) suppressed 
by decoherence due to a momentum uncertainty cTp = 1 GeV. 
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